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Abstract. The computer calculations in [5] to classify sharp polynomials with nonegative 
coefficients constant on the line x + y = 1 have been extended to degree 19, where a 
surprisingly large number of 13 sharp polynomials was found. 



1. Introduction 

Let ~H(2, d) be the space of polynomials p(x, y) of two variables with nonnegative coeffi- 
cients such that p(x, y) = 1 whenever x + y = 1. It is known [1] that the degree d < 2N — 3, 
where iV is the number of positive coefficients of p. Furthermore, it is known that for each 
odd degree d, there exists a group invariant polynomial with precisely d = 2N — 3. We call 
such polynomials the sharp polynomials. See [I], 2} for more information, background, and 
motivation. 

In [2] using computer code, we have have classified the sharp polynomials in 1-1(2, d) up 
to degree 17. Due to the increase in the speed of computers over the last several years and 
improvements to the computer code to be described below, it was possible to run the com- 
puter code for degree 19. The computation with all the improvements takes approximately 
5 days on 4 relatively recent CPU cores. 

The computer has found 13 polynomials polynomials in degree 19 up to symmetry of the 
variables, several of which were unexpected. Two of the polynomials are symmetric. The 
number of sharp polynomials in degree 19 is stunning as there have been only 16 sharp 
polynomials in all odd degrees up to degree 17. 

We have therefore computed a new term for the sequence A143107 on OEIS [3]. That 
is a sequence whose Nth term is the number of sharp polynomials of degree 2N — 3, or in 
other words, the number of polynomials in H(2,2N — 3) with A" nonzero coefficients. In 
this sequence symmetry is not taken into account and therefore there are 24 polynomials 
altogether for degree 19. The sequence is now known to be: 

0,1,1,2,4,2,4,8,4,2,24. (1) 

At the time of this writing, the code is running for degree 21, and it is expected to finish 
in approximately a year. We expect that only the group invariant polynomial is found in 
degree 21, and so we expect the next term in the sequence to be 2. 
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2. The polynomials 
Let us list the polynomials that were found. First we have the group invariant polynomial. 

19 + 19 y x 17 + 152 y 2 x 15 + 665 y 3 x 13 + 1729 y 4 x 11 + 2717 y 5 x 9 + 2508 y 6 x 7 + 

1254 y 7 x b + 285 y 8 x 3 + 19 y 9 x + y 19 (2) 
Then we have several polynomials with integer coefficients. 

+ 19 y x 17 + 152 y 2 x 15 + 665 y 3 x 13 + 1729 y A x 11 + 2090 y 5 x 9 + 627 y 9 x 5 + 

Q27y 5 x 5 + 285 y 8 x 3 + 19 y 9 x + y 19 (3) 

x 19 + 285 y 3 x 13 + 1425 y A x 11 + 19 y 9 x 9 + 2679 y 5 x 9 + 19yx 9 + 2508 y 6 x 7 + 

1254 y 7 x 5 + 285 y 8 x 3 + 19 y 9 x + y 19 (4) 

x 19 + 285 y 3 x 13 + 1425 y 4 x 11 + 19 y 9 x 9 + 2052 y 5 x 9 + 19 y x 9 + 627 y 9 x 5 + 

Q27y 5 x 5 + 285 y 8 x 3 + 19 y 9 x + y 19 (5) 
Next, three polynomials with rational coefficients with denominator 25. 

iq , rt it 2 i. 15371 y 3 x 13 61372/ 4 x n 4807 y 5 x 9 

x 19 + 19y x 17 + 152 y 2 x 15 + + + + 

25 5 5 

1254y 13 x 3 4617 y 8 x 3 1254y 3 x 3 „ iq , . 

5871^ 4617^ 4617^ 

25 5 5 

1254 y 13 x 3 4617 t/x 3 1254 y 3 x 3 in q iq _ 

^ + 5871 ,3. g + 4617j^ + ^ + 4617^ + 
25 5 5 

1254y 13 x 3 4617 y 8 x 3 1254 y 3 x 3 „ iq , . 

Then we have two polynomials with denominator 56. 

iq 855 y x 17 646 y 2 x 15 1938y 3 x 13 2907 y 4 x 11 

x 19 + ^ + y - + 1 + 1 + 

56 7 7 7 

3553 y 5 x 9 323 y 8 x 3 209y 17 x 323 y 9 x 209 yx iq 

19 855 y x 17 323 y 2 x 15 323 y 8 x 9 323 y 5 x 9 

X + 56 + 7 + 7 + 2 + 

323 y 2 z 9 323 y 8 x 3 209 y 17 x 323y 9 x 209 yx iq 

+ + + — + —+y (10) 

7 7 56 28 56 y V ; 
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Finally, rather surprisingly, there are 4 very similar polynomials with denominator 110, 
two of which are symmetric in x and y. 

19 yx 17 323 y 2 x 15 323 y 3 x 13 323 y x 6 



x «. + ^_ + _ f _ + _ f _ + _^ 

,13^4 in. ,17 



323 y Li x A 323y LA x 2 19 y" x 323y b x 399 yx 



+ — - " + — " • " • (11 

II 11 2 55 110 * v 



19 19 yx 17 323t/ 2 x 15 323 y 3 x 13 323 yx 6 
X + ^2~ + H + H + ^5 _ + 



323y 13 x 3 323y 15 x 2 19 y 17 x 323 x 399 „ 

+ + — + — + —+y 12 

11 11 2 55 110 y ' 



1q 19 yx 17 323t/ 2 x 14 323 u 4 x 13 323 wx 6 
19 _| a 1 a 1 y 1 v , 

X + 2 + 11 + 11 + 55 + 



323y 13 x 4 323y 14 x 2 19y 17 x 323 y 6 x 399 yx lq , , 

iq 19 yx 17 323 y 2 x 14 323 y 4 x 13 323 yx 6 
X + 2 + 11 + 11 + 55 + 

323y 13 x 3 323y 15 x 2 19y 17 x 323 y 6 x 399 yx iq . , 

+ + — + — + — + y 9 (14) 

11 11 2 55 110 y y J 

3. The computer code 

The code used is the C code described in detail in [2]. The new revision of the code 
that was used in this computation has been posted at j3]. The basic idea is to consider 
p(x, 1 — x) — 1 =0, which provides a linear equation for the coefficients. We write this 
equation as a matrix A that takes coefficients of degree d—1 or less to the degree d coefficients. 
We know that the degree d coefficients have the form x d + y d . We iterate over the list of 
possible monomials of degree d — 1 or less, taking the corresponding columns of the matrix, 
we look for nonzero solutions. The idea of the algorithm is to find those submatrices that 
are not of full rank (have a nontrivial solution). Then we check if this solution has positive 
coefficients. There are several heuristics that are applied that can avoid doing row reduction 
at all. Already in [2] , to improve speed of the row reduction, we first used mod p arithmetic 
for reduction for a small prime, as most of the submatrices considered are full rank. This 
technique reduces the need to do row reduction in full integer arithmetic in vast majority of 
the cases. We have used p = 19 in degree 17 or less, and in this calculation, we used p = 23. 
The prime must not divide the degree, as most entries in the matrix are divisible by d. 

The major improvements that were done in this revision the following 

(i) Before mode p, the row reduction is first done mod 2. Mod 2 is much faster than mod p. 
The code actually does column reduction in mod 2 and considers columns as unsigned 
integers. The reduction removes one internal loop as adding columns together is simply 
an XOR operation. Unfortunately mod 2 arithmetic only eliminates 90% of the full 
rank cases, but it is approximately 4 times as fast on these 90%. If the test fails on 
mod 2, we move to mod p as before. 
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(ii) Just as the Mathematica code from [2], the C code now ignores polynomials that are 
"right side heavy." That is, polynomials where there are more terms with higher power 
of y than of x. This optimization reduces the run time by approximately one third. 

(hi) Many other more minor optimizations were done, whose individual impact was harder 
to measure. 

Overall, together with improvements in speed of computers since 2008, the current code 
runs approximately 25 times faster than it did in 2008. 
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